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Abstract 

We model the spread of a SIS infection on Small World and random networks using weighted graphs. 
The entry Wij in the weight matrix W holds information about the transmission probability along the 
edge joining node Vi and node Vj. We use the analogy between the spread of a disease on a network and 
a random walk performed on this network to derive a master equation describing the dynamics of the 
process. We find conditions under which an epidemic does not break out and investigate numerically the 
effect of a non-symmetric weight distribution of the initially infected individual on the dynamics of the 
disease spread. 
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1 Introduction 



The study of epidemiological models has been a subject of great interest for 
many years. The aim is to model the spread of a particular infectious dis- 
ease, reproducing the actual dynamics of the disease and designing strategies 
to control and possibly eradicate the infection. Several approaches to tackle this 
problem have been undertaken. The majority of epidemic models are based on 
a compartmental model in which the individuals are grouped according to their 
disease status P, [2]. The basic models describe the number of individuals 
that are susceptible to {S{t)), infected with {I{t)) and recovered from {R{t)) 
a particular disease at time t. The difference in responses between individuals, 
the influence of the topological structure of the system and many other com- 
plex aspects of the progression of the disease are neglected in this approach. 
Although the simplicity of the model means a loss of information and reality, 
it enables us to get a first glimpse of the inner workings of the dynamics of the 
disease spread and makes calculations of threshold values and equilibria pos- 
sible. The assumptions of this model lead to two standard sets of differential 
equations that have provided the foundation of the majority of mathematical 
epidemiology: 



1) The Susceptible-Infectious-Removed (SIR) model: 

^ = i3N-XS-pS 
dt 

^ = XS-dl-pI 
dt 



2) The Susceptible-Infectious-Susceptible (SIS) model: 

^ = 5I-\S 
dt 



(1) 



(2) 



where N is the population size, (3 is the birth rate, p is the natural death rate, 
5 is the recovery rate and A is the infection rate. The SIR model is a suitable 
model for infectious diseases that confer lifelong immunity, for example measles 
or whooping cough. The SIS model is mainly used to model the spread of sex- 
ually transmitted diseases, such as chlamydia or gonorrhoea, where repeated 
infections are common. In these models a random mixing assumption is made: 
each individual has a small and equal chance of coming into contact with any 
other individual. Many modifications to this basic approach have been made to 
account for more heterogeneities. One approach is to further subdivide the pop- 
ulation into subpopulations, with different mixing rates in these groups. This 
means that the parameter j3 in the above equations is replaced by a matrix, 
describing the transmission of infection between different groups. Nevertheless, 



the random mixing assumption, at least within the subgroups, remains un- 
changed. In reahty, however, it is usually the case that the number of contacts 
of an individual is much smaller than the population size and random mixing 
does not occur so that the above model can only serve as a relatively crude 
approximation. Models that incorporate network structure avoid the need to 
rely on the random mixing assumption. They do so by assigning each person 
a fixed amount of contacts. Networks thus capture the permanence of interac- 
tions. A network (or graph) is comprised of a set of nodes and a set of edges. 
An edge is a connection (or bond) that links two nodes. Not all nodes in a 
network are connected directly by one edge. Nodes that are connected by one 
edge are called neighbors. We introduce following network quantities: 

• V = {vi}, the set of nodes. 

• E = {cij}, the set of edges, e^j is the edge running from node Vi to node 

Vj. 

• ki, the degree of node Vi, i.e., the number of neighbors of node Vi. 

• < k >, the average degree of the nodes. 

• kmax, the maximum degree found in the network. 

• P{k), the degree distribution, i.e., P{k) is the percentage of nodes in the 
network that have degree k. 

In a network used to model the spread of a SIS infection, the nodes represent 
individuals that are either infected by or susceptible to the disease under con- 



sideration. Edges represent interactions between individuals: The disease can 
only spread from one node to the next if there is an edge connecting them. If 
we assign a weight Qij^ representing the probability of node Vi to infect node Vj^ 
to each edge e^j, we obtain a weighted graph. The weight matrix W = {wij] 
is defined as follows: 



This matrix gives information about the network topology, the relations be- 
tween individuals in it and the characteristic of the disease under consideration. 
The use of networks in Mathematical Epidemiology has grown exponentially 
since the middle of the 20th century. In recent times, starting with the work 
of Pastor-Satorras and Vespignani , there has been a burst of activity on un- 
derstanding the effects of the network topology on the rate and pattern of the 
disease spread [[3]-[7]]. Amongst others the main network types studied are the 
following: 

1) Regular lattices 

In regular lattices each vertex is connected to its k nearest neighbors, to form 
either rings (one-dimensional) or grids (two-dimensional) [8]. 

2) Random graphs 

The term random graph refers to the disordered nature of the arrangements 




Qij^ if node Vi is a neighbor of node vj 







else. 



of links between different nodes . Erdos and Reyni (ER), in their first paper, 
proposed a model to generate random graphs with N nodes and K links [9]. 
Starting with N disconnect nodes, these random graphs are generated by con- 
necting couples of randomly selected nodes, prohibiting multiple connections, 
until the number of edges is K. 

An alternative model for ER graphs is created by connecting each pair of nodes 
with probability < p < 1. This procedure results in graphs having differ- 
ent amount of edges present but the two models show strong similarities and 
coincide in the limit of large N. random graphs have degree distribution ap- 
proximately Poisson with parameter < k >. 

3) Small World networks 

The study of dynamical processes over real networks has pointed out the ex- 
istence of shortcuts, i.e. links that connect different areas of the graph, thus 
speeding up the communication between otherwise distant nodes. This is known 
as the Small World property . It is mathematically characterized by a relatively 
short average path length that depends at most logarithmically on the network 
size. This property is observed in a variety of real networks including random 
graphs. To distinguish between random and Small World networks, the Small 
World property is often associated with the presence of clustering. Watts and 
Strogatz have proposed to define Small World networks as networks having 
both, a short average path length (like random graphs) and a high clustering 



coefficient [TO]. 

4) Scale free networks 

Networlcs tliat liave a power law degree distribution: P{k) ~ k~^, 2 < 7 < 3 
are called scale free [TT]. In comparison, in a random graph, P{k) decays faster 
than exponentially. In scale free networks we often find a significant amount 
of nodes with very high degree. In the context of disease spread on a graph, 
these nodes are aptly called super spreaders. In scale free networks the average 
degree < /c > is no longer a relevant variable and one expects fluctuations in 
< > to play an important role. 

In their work, Shirley and Rushton investigated how the speed of the disease 
spread is influenced by certain topological characteristics of the graph [12]. 
They found that an epidemic spreads fastest on a scale free network, followed 
by random graphs and is slowest on regular lattices. Small World graphs lay 
in between random and regular graphs. In their model, the transmission prob- 
ability between nodes is homogeneous throughout the network, i.e., wij = w. 
However, when modelling the spread of a disease on a network one should, along 
with complex topological features, take into account heterogeneity in the inten- 
sity strength between nodes: not every individual is susceptible to infection 
or capable to infect its neighbors to the same degree. The complexity in the 
capacity and intensity of the connections also plays an important role in other 
real networks like scientific collaboration networks, air-transportation networks. 



internet clusters and other large infrastructure systems. Newman showed that 
weighted networks can in many cases be analyzed by using a simple mapping 
from the weighted graph to an unweighted multi graph and then applying stan- 
dard techniques for unweighted networks [13J. Simonsen et al. studied diffusion 
on unweighted and weighted networks and discovered that the eigenvalues of 
the transfer matrix describing the process can be used to recover large scale 
topological features of the system [I4]-[15]. Most recently, Vasquez introduced 
a type-network representation of a Small World weighted graph to take into 
account the population heterogeneity in a very general approach [16]. He ob- 
tained a recursive equation for the probability distribution of the outbreak size 
as a function of time and demonstrated that the expected outbreak size and 
its progression in time are determined by the largest eigenvalue of a certain 
matrix (reproductive number matrix) and the characteristic distance between 
individuals. 

In this work, we focus on the effect of heterogeneities in the transmission prob- 
abilities on the dynamics of the disease spread. We use the strong analogy 
between the spread of a disease on the network and a random walk performed 
on that network to derive a master equation describing the dynamics of the 
process. We find conditions under which an epidemic does not break out and 
investigate numerically the effect of a non-symmetric weight distribution on the 
dynamics of the disease spread. 

This paper is organized as follows: In section 2, the epidemic models (SIS) is 



described. Section 3 contains the derivation of the master equation describing 
the dynamic process. Section 4 is concerned with identifying conditions that 
prevent the outbreak of an epidemic. In section 5, we discuss our numerical 
simulations and their results. We close with an overview of the main results 
and a discussion in section 6. 

2 The model 

In this paper we will consider the standard SIS epidemic models without birth 
and death on a Small World graph. 

2.1 The epidemic model 

At each time t, the population is divided into two categories: susceptibles, 
Sit) > 0, and infectious, I{t) > 0. We normalize so that S{t) + I{t) = 1 
for all times t. Susceptible members are virgin territory for the disease, whilst 
infectious members are both infected and capable of infecting others with whom 
they are in direct contact, i.e., their neighbors. After being infected for time 
Tinf, an individual returns to the susceptible class. In this paper we take 
rinf = 1- 

2.2 The Network 

In our networks, the vertices of the graph represent the individuals of the pop- 
ulation under consideration, the edges describe the contact patterns between 
individuals. The adjacency matrix A of this network is defined by 



j 1, if node Vi and vj are joined by an edge 
aij = < 

I 0, else 
and gives information about neighbor relations in the network. Assigning a 
weight Wij to the edge connecting node Vi and node Vj, representing the proba- 
bilities for the disease to spread from node Vi to node Vj, will result in a weighted 
graph giving information about 

1. The network topology describing the population and the relations between 
individuals in it and 

2. The characteristics of the disease under consideration. 

This information is encoded in the weight matrix W, defined in the previous 
section. 

The weights play the important role of conveying susceptibility and transmis- 
sibility levels of individuals. For example, consider nodes Vi and Vj with Wij 
being relatively large in comparison to wji. This could either mean that node 
v'^s capability to transmit the disease is very large and VjS very low or that node 
Vi is not very susceptible to transmission of the disease from node Vj while Vj 
can easily catch the disease from Vi. 



3 Random Walks and Epidemics 

In this section we will describe the model implemented to simulate the spread 
of an infectious disease throughout a population and derive a master equation 
governing the dynamics of the system. The design of the model was inspired by 
the work of Alves et al [17]. There is an strong analogy between the model of 
an epidemic on a network and a random walk performed on this network that 
can be described as follows: 

Suppose we want to follow the spread of an SIS epidemic on a social network. 
The individuals of the network are represented as nodes. At each node we 
place a certain amount of random walkers. Each random walker represents 
the possibility of an infection to happen: If a walker moves along the edge Cij 
from the infected node Vi to the susceptible node Vj, node vj has been infected 
with the disease. If no walker moves to vj, the node remains susceptible. The 
probability of the walker moving from Vi to vj is given by the entry Wij of the 
weight matrix W. There is an artificial component to this model: The number of 
walkers placed at each node at time to depends on the length of our experiment 
and the infection time Tinf. Since each walker represents the possibility of one 
infection, we need to have a sufficient amount of walkers at each node at time 
to so that we do not 'run out' of walkers (the ability to infect) at some time 
t before the end of our experiment. If we want to model the spread of a SIS 
infection for T^ax time steps, we need to place at least walkers at 

node Vj. This is to ensure that, when the node is infected, an infection can 



happen along each bond eminating from our node throughout the duration of 
the experiment. 

Let us define following quantities: 



(Here, the sum is taken over all nodes.) 

• Let W = {wij} be the matrix of transmission probabilities described ear- 
lier: Wij is the probability of infected individual Vi to infect susceptible 
individual Vj. 

• Let l{t) = {liit)}, where 



I{t) is called the infection matrix. The number of non-zero elements of I{t) 
is the total number of infected individuals at time t. 

We model the spread of the infection through the population via a random walk 
performed on the network: 

At t = to, we place Tmax^i walkers at each node Vi. Each walker represents 
the possibility of an infection to happen. Choosing T^axki walkers (and hence 
= Tmax Ylj kj) ensures that no node looses the ability to infect its neighbors 
(provided the node is infected) before the maximum number of time steps T^ax 
(i.e., the end of the experiment) is reached. At each time step t, at a generic 
node Vi, one of two things may happen: 



Let N = Tmax X^j kj be the total number of random walkers participating. 




1. If node Vi is infected: ki walkers may each make one move: The walkers 
are allowed to move between adjacent vertices. What edge, out of the 
possible outgoing ones, a walker chooses to move along is picked at random 
with probability equal to the weight assigned to this directed edge. If the 
walker moves to a susceptible neighboring node, that node will be infected 
with the disease. Walkers from infected neighboring nodes may enter. After 
one time step the status of Vi is reset to susceptible (/^(t — 1) = 1, Ii{t) = 0). 

2. If node Vi is susceptible: No walkers leave Vi but walkers from neighbor- 
ing infected nodes may enter. As soon as a walker enters the susceptible 
node Vi, it is considered infected {Ii{t — 1) = 0, = 1). At the next 
time step ki walkers of node Vi may move to infect any of the susceptible 
neighbors of node Vi. 

Let us denote by r]i{t) the percentage of walkers at node Vi at time t. This is 
the quantity that will give us information about the dynamics of the disease 
spread. We note following properties of r]i{t): 

• If node Vi is susceptible at time t, r]i{t+l) —T]i{t) >= and the larger this 
difference is, the more susceptible to infection is the node. 

• If node Vi is infected at time t, walkers are leaving but may also enter from 
neighboring nodes so rii{t+ 1) — rjiit) does not give any useful information. 

• The change in walker density from time t to time t + 1 for the susceptible 



node Vi satisfies: 

kr 



< hit + 1) - r]i{t)\ < J2 



k- 

If 

0<|ry,(t + l)-ry,(t)|= 

mGNeighi ^3 

all neighboring nodes are infected and all walkers from those nodes move 
to Vi. 

• The change in walker density from time i to time t + 1 for the infected node 



Vi satisfies: 



If 



all ki walkers leave node Vi and no walkers enter. 



• Hence, for a generic node i'i, we have 



o<te(* + i)-'?K*)l < ^ ^ + *=■ 



max 



m&Neighi ^3 ^3 

T~[^max ~l~ ^max] (3) 

• Nodes that are extremely susceptible to infection at time t have large \r]i{t+ 
l)—'r]i{t) \ values (>> y^V). Infected nodes that satisfy \r}i{t+l)—r}i{t)\ ^ 
j^'k 1 ^^^y infectuous. Hence \rji{t+ 1) — rii{t) \ quantifies the extend to 
which node Vi participates in the spread of the disease under consideration. 



Since the infection time Tin/ is one and \r]i{t + 1) — r]i{t)\ < ^jj^^j^. + 



for Tjnax > 3, we have 



1 



I f-> I J'^max I I r, I y'^max J 



and for T^.:. > 1: 



T T 

y-'- max I £-< iJ'^max ' \_ c 



— Vii'-^max) 



• From the derivation of these bounds, we can conclude that a generic node 
Vi satisfies: 
1 



2 2 J ^^^kr^ 



• Nodes with large r]i{Tmax) values have been infected with great intensity 
but have transmitted the disease with low frequency. 

• Nodes with small r]i{Tmax) values have rarely been infected themselves but 
if infected have contracted the disease to many of their neighbors. 



E,r/.(t) = lfor alH. 



3.1 The Master Equation 

We now derive the equation governing the dynamics of the disease spread. This 
equation is called the master equation. The change in the walker density is the 
difference between the relative number of walkers entering, J^{t), and leaving, 
Jj^(t), the same vertex over the time interval t — > t + 1. Hence, 

^See appendix A and B 



ij,{t + l) = ij,{t) + Jr{t)-J+{t). (4) 

For the moment, let us presume that node Vi is infected. Then the edge current 
on the directed edge from vertex Vi to a neighboring vertex vj is given by: 



im 



(5) 

The edge current is the fraction of walkers moving along this edge according 
to the weight distribution emanating from that node. Note that ^ 



is the probability of a walker deciding on the edge from vertex Vi to vertex vj. 
Since no walkers are leaving if the node is susceptible, we must have: 



jrit)= (6) 

and 



jQNeighi 

and upon substitution into equation (jlj), we obtain 



(7) 



j^Neighi j€Neighi 
j^Neighi 

With the matrix T defined as follow 



J- max Z^meNeighj '^jm ^ max 



if Vj is a neighbor of vi 

else, 
we can rewrite equation (EI) in matrix notation. 



J- max J- max 

t t 



= r,(t,) + ^(Y^ I(r)) . T„(«„) - I(r)) • „(«„) 

J- max ^ , ^ , 

T=to T=to 



T 

-'- max , 

T^to 



Equation ([9]) is the desired master equation. This equation governs the dynam- 
ics of the disease spread. We can see that r]{t+ 1) depends on the initial walker 
distribution, the sum of the infection matrices from time to t, which tells us 
which nodes have not been infected yet and the matrix T, which is called the 
transfer matrix. 



4 Conditions for non-outbreak 



If 'r}{t + 1) « ^(^o) for all t, an outbreak of the epidemic does not take place. 
From equation (j9]) we can see that this is the case if either 

or 

In this paper we will concentrate on the first condition. We investigate how 
the distribution of incoming and outgoing weights of the initial node effects 
the dynamics of the disease spread. We define a weighted difference Di and 
a nodal entropy Si to quantify the difference between incoming and outgoing 
transmission probabilities of a generic node Vi. There are two ways in which 
condition 1 can be realized: 

1) The epidemic does not break out: either non or relatively few neighbors of 
the initially infected node get infected by it and they, in turn, do not infect a 
significant number of their neighbors until, after a short number of time steps, 
there are no infected individuals in the population. This happens if the majority 
of the transmission probabilities stay below the threshold value of the network. 

2) Infection of susceptible individuals takes place over a long period of time but 
infected individuals stay localized around the initially infected node. 

To investigate when this may happen, we introduce the following difference 
measure: For node Vi, we define 



i.e., Di describes the difference between outgoing (wij) and incoming {wji) 
transmission probabilities of node Vi. If Di < 0, the sum of the outgoing 
weights is smaller than the sum of incoming weights and vice versa if > 0. 
Nodes with Di < are more likely to be infected often but infect their neighbors 
with low intensity. Nodes with Di > are more infectious but less susceptible 
to infection. We expect some sort of localization of the disease spread around 
the initial node Vi to occur if Di <C and the average of the weights of the 
other nodes are sufficiently small. Let us also define an entropy describing the 
diversity of the differences of incoming and outgoing weights of a node: 

Let 

Pi{d) = % of edges eij of node Vi with Wij — Wji = d 

Then Pi{d) is a probability distribution on the differences of incoming and 
outgoing weights of node Vi (i.e., '^(iPi{d) = 1 for all i) and we utilize this 
distribution by defining the following nodal entropy in the standard way: At 
node Vi, we let 

S^ = -J2PiWo9iPi{d)). 

d 



gives information about the disorder of these differences: The larger the 
value of Sj^, the large the variety of differences in incoming and outgoing 
weights. 

5 Numerical investigation on a Small World Graph 
5.1 The network 

We create a Small World graph in the following way: N points lying in the 
square [0, 1] x [0, 1] are randomly selected. These points represent our popula- 
tion. We choose the following network defining quantities: 

• r: The short distance radius 

• pr'. The probability of short distance bond formation 

• R: The long distance radius 

• Pr: The probability of long distance bond formation 

• Wr'. The transmission probability along short distance bonds. 

• Wr: The transmission probability along short distance bonds. 

Points lying a distance less than r away from each other are connected with 
probability Pr- Points lying a distance more than R away from each other are 
connected with probability pij. Transmission of the disease along short distance 
bonds occurs with probability j?^-, transmission along long distance bonds occurs 
with probability Pr. 



5.2 Sensitivity to initial conditions 

To see when Yll-=t^ I(t) ~ 0, let us set up our Small World graph as follows: 
The population size is 2000. All transmission probabilities across bonds, 
regardless of the nature of the bond, except for node one, are held fixed at 
3% (i.e., Pr = Pr). This value is chosen because our numerical simulations 
with homogeneous transmission probabilities showed that an epidemic does 
rarely break out for values around 3% We will look at two different scenarios 
to evaluate the influence of long distance connections in the network: 

1. Pr is chosen small enough so that the average long distance degree of a 
node Hes at 0.09 (Small world graph). 

2- Pr = 0, no long distance connections are permitted (random graph). 

Pr is chosen so that the average short distance degree is 25. Node one, the 
initially infected, has degree ki = 20 with 19 short distance connections and 
one long distance connection in case one and 20 connections in case 2. We deflne 
the participation ratio PR{t) at time t as the percentage of the population that 
has been infected at least once during the time interval [0, t] and the diameter 
Diam{t) as the maximum distance from the initially infected node where an 
infected individual can be found during the time interval [0,t]. We investigate 
how the diameter and the participation ratio change with varying Di and Si 
values, i.e. how the non-symmetry of the transmission probabilities of the initial 
node and the extend to which it is not symmetric effect the dynamics of the 
epidemic. The diameter and participation ratio are measured at the last time 



step Tmax- Tmax is chosen large enough so that either the epidemic has spread 
through all the population, resulting in a maximal diameter and participation 
ratio, an equilibrium is reached (participation ratio and diameter approach 
some limit) or the epidemic has died out for some t < T^ax- Our numerical 
simulations resulted in figure number 1. 

We can see that the plots look very similar for graphs without long distance 
connections and Small World graphs. As expected, the introduction of long 
distance connections between nodes drives up the diameter and participation 
ratio. 

We note that the plots of the diameter and participation ratio in both cases 
are almost symmetric about the line Di = 0. For small values of 5f we 
observe a basin around Z)i = 0: diameter and participation ratio stay low 
as we hoped. But note that, as 5f increases, the behavior around Di 'P^i Q 
becomes increasingly erratic and unpredictable. As we move away from Di ~ 
regions but stay in areas where Si is relatively low, the diameter increases. 
We can see that the infection reaches the outskirts of our domain [0, 1] x [0, 1]. 
participation ratios climb as high as 40%. Again, as S'f increases we notice the 
smoothness of the plots breaking and see disorder taking over. Looking at the 
participation ratio plots we notice an oddity: At Di « 0.4 and Si ~ 0.65 we 
see a sharp spike much higher than any surrounding participation ratio values. 
This behavior is surprising and needs further investigation. 



6 Results and Discussion 



We have used the analogy between the spread of a disease on a network and a 
random walk performed on this network to derive a master equation describing 
the dynamics of the process. We found two conditions under which an epi- 
demic does not break out. One of these conditions strongly depends on the 
transmission probabilities of the initial node. This lead to the consideration 
of a non-symmetric weight matrix W. The majority of research concerning 
epidemic modeling on networks assumes that the transmission probability ma- 
trix W is symmetric. This means that node Vi infects neighboring node vj 
with the same probability as node vj infects node Vi. We have focused on a 
more realistic setting where we take into account the heterogeneity in the in- 
tensity strength between nodes, i.e., a non-symmetric weight matrix W. In 
particular, we focus on the initially infected individual showing these hetero- 
geneities, leaving all other nodes with homogeneous transmission probabilities. 
We chose these small enough so that in case of a fully homogeneous network 
an epidemic would not break out. We numerically investigated the effect of 
the non-homogeneous weight distribution of the initial node. To quantify the 
heterogeneity in transmission probabilities, we defined two nodal quantities: 

1. The difference Di = ^^^^^'^^jj — 'A^ describes the difference between out- 
going {wij) and incoming {wji) transmission probabilities of node Vi. 

2. The entropy Sf^ = — ^dPi{d)log{Pi{d)), where 

Piid) = % of edges of node Vi with Wij — wji = d, 



describing the disorder in the differences d. 

Through our numerical simulations we obtained surface profiles of diameter 
and participation ratio for a variety of Di and Si values for graphs with long 
distance connections and without. Graphs without long distance connections 
are just random graphs. Introduction of long distance connections creates Small 
World graphs. We obtained the following results: 

• The quasi-symmetry around Di ^ for low Sf^ values shows that regardless 
of whether incoming weights of node 1 are larger than outgoing or vice versa, 
the dynamics of the disease spread is very similar. 

• An increase in <S*f values leads to a breakdown of symmetry and unpre- 
dictable dynamics. This breakdown of symmetry occurs at values about 
0.5. 

In summary, we have shown that the introduction of even very few non-symmetric 
transmission probabilities along edges of an otherwise homogeneous network 
changes the dynamic of the disease spread. For a larger disorder in the dif- 
ferences between incoming and outgoing weights of the initially infected, the 
dynamics becomes unpredictable. The threshold value of this disorder measured 
by the entropy 5f lies at about 0.5. It is quite apparent that a realistic model 
of the spread of an infectious disease throughout a population should take into 
consideration the heterogeneities in transmission probabilities between individ- 
uals. We have taken a step towards this direction by investigating the effect of 
the introduction of a single node with heterogeneous transmission probabilities. 



Further research opportunities present themselves. One should investigate the 
effect of heterogeneities in more than one node. One should try and find a mea- 
sure quantifying the extend of anti-symmetry in the network as a whole and 
investigate if there exist certain threshold values determining if an epidemic 
breaks out or not. We also plan to analyze the second condition for a non- 
outbreak found in this paper. This will involve a thorough investigation of the 
properties of the transfer matrix T. Another direction of research would be to 
run the simulation on scale-free networks to see if we find the same symmetry 
about Di = 0, if a breakdown of this symmetry occurs and if so, what the 5f 
threshold value is. 

Network analysis in the context of epidemic modelling has helped us create more 
realistic settings to investigate the spread of infectious diseases throughout a 
population. The study of non-symmetric weight matrices may prove useful 
in creating even more realistic models taking into account the variability of 
transmission probability between individuals of the population and could lead 
to more accurate threshold values predicting the outbreak of an epidemic, hence 
enabling us to develop better methods to prevent or eradicate the disease at 
hand. 
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7 Appendix 
7.1 A 

Derivation of lower bound 



y-'-max \ ^ IJi^maxi L o J — 'hy-'-max) 



First, we make some observations: 

Tinf = 1 means that a generic node Vi can be infected by its neighbor at most 
every second time step. Since a node can only loose walkers when it is infected 
and node Vin is the only infected individual at time ^o, wc conclude that Vin is 
the node that potentially has the least number of walkers at time Tmax- We 
also need to assume that node Vm is of maximum degree kmax- The following 
chain of events will result in node Vin having the smallest number of walkers at 
the end of the experiment: 

• Time to: Vm is infected, Vj is susceptible. 



/ N m 



max 

max'' 



l^m ^' 

• Time ti: All km = k^ax walkers leave Vm- 

k 

(tl) — {Tmax — 1)- 



m 



max 



t'inturns susceptible. 

• Time 12'. To keep loosing walkers at the largest rate, the node needs to 
be infected at every second time step. We assume that only one walker 



causes the infection at each time, keeping the amount of walkers entering 
at a minimum. 

^m(^2) ~ {Tmax 1) t^^T 1" 

• When infected, the node looses all kin walkers in the next time interval. 

• When susceptible, one walker enters, causing infection. 

• Hence, for t = Tmax, we obtain: 

n- (T ) - ^ \(T r ^rrrnx ny , 7^ , 

'liny-'- max) — i. y-'-max \ 2 IJ'^max ^ L 2 



7.2 B 

Derivation of upper bound 



r^.(T ) < ^ (T r (^max 3) r jTmax 2) 3' 

iliy-'-max) _ , \-i-max \ l^'^maa;^^ | | V'^max^ 



Let be a neighbor of the initially infected individual Vin- If we assume that Vj 
and all its neighbors are of maximum degree kmax and notice that only walkers 
sitting at infected nodes can move, we conclude that the most walkers can 
accumulate at a node Vj, neighbor to vin. The following chain of events will 
result in Vj obtaining the largest number of walkers over the time span of the 
experiment: 



Time to: Vm is infected, Vj is susceptible 



max'' 



^m 



l^m '^r. 

Time ti: All walkers leave node Vin and move to node Vj. Node Vj gets 
infected. Node Vin turns susceptible. 

rijitl) = iTmax + l)^^. 

l^m ^ 

Time i^'. All k^nax walkers leave node vy. 

h 

r]j[t2) = ylmax)- 



m 



infecting each neighbor. The node returns to the state of susceptibility . 

Time ^3: Assuming that all walkers of all infected neighboring nodes will 
move to Vj, we obtain: 

^m '^rn .^^^ l^rn '^rn 

/rp \ ^max {k-max) 

— [J-max)^ h \^ h ' 

since we assumed all neighbors of Vj to be of maximum degree. 

We now assume that step 2 and 3 repeat until the end of the experiment. 

Time t^: 

/rp -1 \ ^max . ij^maoc) 

V-^ max J- j ] r 



m 



Time ^5: 

/rp -1 \ ^max . r> ij^max) 

\^ max -L j r ^; 



m 



K-'-max I D U'^maa; ~r | 



I'm 
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